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Abstract 

We present a precision determination of the critical coupling fi c for the deconfinement transition in pure SU(2) gauge theory in 2 + 1 
f-*) dimensions. This is possible from universality, by intersecting the center vortex free energy as a function of the lattice coupling 
CXl j6 with the exactly known value of the interface free energy in the 2D Ising model at criticality. Results for lattices with different 

numbers of sites Af, along the Euclidean time direction are used to determine how ft varies with temperature for a given N, around 

the deconfinement transition. 

'Key words: Center vortex free energy, twisted boundary conditions, deconfinement transition, universality, 
j-^, PACS: 12.38.Gc, 12.38.Aw, 11.15.Ha. 
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TT"1- Introduction 

a : 

I It is widely accepted today that the deconfinement transition 
Onin pure SU(A0 gauge theories at finite temperature is driven by 
^ the dynamics of center vortices HI . In the vortex picture of con- 
i finement, Wilson loops acquire a disordering phase factor from 
every vortex that they link with. The area law for timelike Wil- 
son loops in the pure gauge theory comes from the percolation 
of spacelike vortex sheets in the confined phase. Their free en- 
ergies have been measured over the deconfinement phase tran- 
sition at finite temperature in the 3 + 1 dimensional pure SU(2) 
gauge theory from ratios of partition functions with 't Hooft's 
twisted boundary conditions in temporal planes, forcing odd 
numbers of Z2 center vortices through those planes, over the pe- 
riodic ensemble with even numbers A Kramers- Wannier 

duality is then observed by comparing the behaviour of these 
center vortices with that oft Hooft's electric fluxes which yield 
the free energies of static charges in a well-defined (UV-regular) 
way [4], with boundary conditions to mimic the presence of 
'mirror' (anti)charges in neighbouring volumes. This duality 
follows that between the Wilson loops of the 3-dimensional Z2- 
gauge theory and the 3D-Ising spins, reflecting the universality 
of the center symmetry breaking transition. Here we study the 
vortex free energies of pure SU(2) in 2 + 1 dimensions over 
the deconfinement transition because the relevant interface free 
energies of its universal partner, the self-dual Ising model in 2 
dimensions, are all known analytically. Moreover, the vortex 
free energies in 2 + 1 dimensions are much cheaper to simulate 
and discretisation effects vanish more rapidly than in 3 + 1 di- 
mensions. Together these reasons allow for numerical studies 
of much higher precision. 

In Sections 2 and 3 we briefly outline the basic concepts and 
our numerical procedure. Our results are presented in Sec. 4. 
These include the precise determination of the critical coupling 
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for various lattices with up to N, = 9 sites in the Euclidean 
time direction, an analysis of the finite volume corrections and a 
brief comparison with the corrections to scaling in the 2D Ising 
model. We then determine how the critical coupling depends 
on N,, including 1 /N t corrections, and use this result to derive 
how the lattice coupling /? varies with temperature for a given 
fixed N, around the deconfinement transition. This is needed, 
for example, for a detailed finite-size scaling analysis of the 
vortex free energies in 2+1 dimensional SU(2) |5|]. 



2. Concepts and Methods 

For pure SU(2) gauge theory, 't Hooft's twisted boundary 
conditions fix the total number of Z2 vortices modulo 2 through 
each plane of a finite box [6]. Twist in a plane corresponds to an 
ensemble with an odd number of Z2 vortices through that plane. 

In a L 2 x l/T Euclidean box, we can distinguish between 
two types of twist. Magnetic twist is defined in the purely spa- 
tial plane and forces vortices that run along the temporal direc- 
tion. They may spread independently of the temperature T and 
play no role in the deconfinement transition: their free energy 
is expected to vanish in the thermodynamic limit at all temper- 
atures which has been demonstrated explicitly in 3 + 1 dimen- 
sions 1 7]. On the other hand, vortices from twist in the two 
temporal planes are squeezed more and more as T is increased. 
They may no longer spread arbitrarily and this is what drives 
the phase transition. 

In this paper we're interested in configurations with a twist 
in one of the temporal planes. If we denote the partition func- 
tion of this ensemble by Z tw (L, T), its free energy per T is de- 
fined via the ratio 

Z tK (L,T)/Zo(L,T) = e- F >^ (1) 
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Figure 1: Universality 




where Zq is the partition function of the periodic ensemble. 

As we approach the thermodynamic limit, Z m ,/Zo converges 
to a non-trivial fixed point at the critical temperature. It can 
therefore be used as a phenomenological coupling in the same 
manner as the Binder cumulant. Typically, one uses the pair- 
wise intersections of Binder cumulant curves from different lat- 
tice sizes to estimate the critical coupling or temperature [8, 9]. 
Using a ratio of partition functions in this way was proposed by 
Hasenbusch for the 3D Ising model ifioll . 

Here we can go a step further. Since SU(2) in 2+1 dimen- 
sions is in the same universality class as the 2D Ising model, the 
the value of Z,„ ./Zq at the fixed point is exactly known. A single 
temporal twist in SU(2) corresponds to a single anti-periodic 
direction on a N x N Ising lattice. This forces an odd number 
of spin interfaces perpendicular to the anti periodic direction. 
The corresponding ratio of partition functions Z ap /Z pp in the 
Ising model gives the free energy of the system in the same 
way as equation ([T]). In both cases, configurations of minimum 
action/energy dominate the partition function. Thus, in the ther- 
modynamic limit, the free energy of a single spacelike vortex in 
SU(2) is identified with the free energy of a single interface in 
the 2D square Ising model at the respective critical points. This 
universal value is given by lllltl 



lim Z ap (T c )/Z pp (T c ) = 1/(1 +2 3/4 ). 



(2) 



For SU(2) on the lattice, l/T — N,a where N, is the number 
of sites in the time direction and the lattice spacing a = a(J3) 
depends on the coupling. So the critical temperature T c corre- 
sponds to a critical lattice coupling f5 c oo for each N,. The sub- 
script oo reminds us that we only have a strict phase transition 
and hence critical coupling in the limit of infinite spatial vol- 
ume. Still, the intersection of Z, w /Zq with the universal value 
(O gives a reliable estimate of p c ,oo{Ni) provided that the spatial 
length L = N s a is large enough. This will be more precise in 
general than the estimates obtained via pairwise intersections Q 

We assume a finite size scaling (FSS) behaviour of the form 



Z fw /Z = 1/(1 + 2 3/4 ) + b{J3-p c )N] lv + cN; w + . 



(3) 



where u is a correction to scaling exponent that should be ap- 
proximately independent of A^,, and v = 1 is exactly known 
from the 2D Ising model. Furthermore, fi c = /3 c ,oo(N t ), and we 
define a 'pseudo-critical coupling' in a finite volume, /3 c (N t , N s ), 
by the requirement that the corrections to the universal value in 





1.04 - 




1.02 - 






a 


1 - 






T 




0.98 - 


a 


0.96 : 




0.94 - 




0.92 - 



4 universal value 
N„ = 20 
N s = 48 
N s = 96 



8.02 8.03 8.04 8.05 8.06 8.07 8.08 8.09 8.1 



Figure 2: Fitted curves for N t = 5, selected for clarity. 

0} vanish. These estimates then converge to the infinite volume 
critical coupling /J c> oo(AQ as 



/3 c (N t ,N s ) = PcANd-diNdK 



-(oM-l/v) 



(4) 



Here, the coeffiecient d — c/b is an N t dependent fit parameter. 

For notational simplicity we will hereafter drop the sub- 
script oo on the critical coupling. Unless N s is explicitly men- 
tioned, f3 c refers to the infinite volume limit. We have used a 
variety of spatial lattice sizes and the above scaling formula (0) 
to first find p c {N„N s ) for N, = 4, 5, 6, 7, 8 and 9 with high 
precision, in order to then determine the corresponding f3 c (N t ) 
from fitting the volume dependence of /3 C (N,, N s ) by 



3. Numerical recipe 

Twist was implemented in the usual way lfl3ll . We kept pe- 
riodic boundary conditions but flipped the coupling f3 — > -J3 of 
a stack of plaquettes perpendicular to the plane of the twist. In 
other words, we introduced a Z2 Dirac string that is closed by 
lattice periodicity. The result is a transformation of the usual 
Wilson action. 

In practice, the overlap of Z, w and Zq is poor. To overcome 
this we interpolated in the number of flipped plaquettes using 
the snake algorithm of Ref . Ill4ll . This was combined with the 
other variance reduction tricks therein. 

For each combination of N, and N s we performed simula- 
tions for ~ 10 values of /3 around the intersection of Z m /Zq 
with the exact Ising value (fJJ. Random errors were estimated 
via the boostrap method. Since SU(2) lattice gauge theory is 
less computationally expensive in 2+1 dimensions than in 3+1 
dimensions, we were able to perform a very large number of 
measurements. 1-30 million configurations were used for each 
y8, depending on the lattice size. 

It turns out that the free energy F, H ,(J3) = -lnZ, M ,/Zo has 
less curvature than Z, M ,/Zo near the critical point, so it's a better 
candidate for linear approximation. Therefore we performed 
least squares fits, with parameters f\ and /3 C , of the form 



'We thank M. Hasenbusch for pointing out Ref. 1 2|. where a similar idea 
was applied to the Z2 gauge theory in 2+1 dimensions. 



^ w ( J 0)=/i(y8-y6 c ) + ln(l+2 3/4 ) 



(5) 
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Figure 3: FSS data collapse for 16 2 X 4, 24 2 X 4 and 32 2 X 4. 




Figure 4: All values for p c (N t , N s ) obtained in this work. The error bars are 
generally smaller than the symbols, and the surface was rendered from the fits 
to 0} with d(N t ) = 0. 134 N? +1 - 0.5 1 N? according to {§) with d t = 0. 



to obtain estimates fi c = j3 c (Nt,N s ) for the critical lattice cou- 
pling. In each case the reduced ;f 2 was ~ 1, which justifies the 
linear ansatz (O. See Fig. [2]for some representative fits. 

Note, however, that a small reduced^- 2 does not exclude the 
existence of significant systematic errors in /3 c (Nt,N s ). On fi- 
nite lattices, F, W (J3) has positive curvature near the critical cou- 
pling, so p c tends to be underestimated. To control this, we 
carefully chose the size of our fitting windows. For each A^ f , we 
performed precise measurements of F M in a quadratic fitting 
window for one or more of our smallest lattices. By writing F, w 
as a function of the finite size scaling variable 



x = N s t oc +L/g ± 



(6) 



where t = T/T c — 1 is the reduced temperature and £ — ^2|fr v 
are the correlation lengths for T <; T c with v = 1 for the 2D 
Ising model, we were able to translate this data to larger N s . 
See Fig. [3]for an example of FSS data collapse in a large win- 
dow for several lattice sizes@ The relationship between tem- 
perature and coupling is described in Section 14.41 It requires 
a paramerisation of f3 c vs N t - l/T c a c , which we roughly ob- 
tained using literature values of the critical coupling and our 
own preliminary results. The translated data was used to es- 
timate the slope and curvature of F tw (J3) near the critical point 
for large lattices without performing additional simulations. We 
then adjusted the linear fitting windows such that the systematic 
errors in f3 c {N,, N s ) should be less than one quarter of the quoted 
random error for each extrapolated f3 c (N,). 



4. Results 

4.1. Determining p c (N t ) 

We have obtained estimates of /3 c (N t ) for every N, between 
4 and 9. In each case we used 8 spatial lattice sizes, keeping 
an aspect ratio of approximately 3:1 for the smallest lattice and 
using a maximum lattice size of 96 2 x N,. See Fig. |5]for plots 



of fi c (N,,N s ) vs N s for N, = 4, 5 and 6. The data was fitted ac- 
cording to the equation with all three parameters free. The 
plots for N, =1,8 and 9 look very similar. Note the rapid con- 
vergence to the infinite volume values, which Hasenbusch also 
observed for the pairwise intersection of Z ap /Z pp curves in the 
3D Ising model 1 1 Ofl - He found much more rapid convergence 
than for the intersections of Binder cumulants. 

We summarize our results in Table [T] For each N, we in- 
clude the estimates /3 C (N,,N S ) from our two largest lattices as 
well as the fitted values of /3 C (N,). It's clear from the reduced 
^ 2 s that the data is very well described by the FSS ansatz ©. 

For N t = 4, 5 and 6 we were able to surpass the precision of 
current literature values of the critical couplings by two orders 
of magnitude. These lattices also gave us good precision for 
the correction to scaling exponent a>. Our results for N t = 7, 
8, 9 are somewhat less precise, especially the fitted values for 
to. This is because the lattices used had smaller aspect ratios 
and we collected fewer statistics. Still, to is consistent between 
each of the fits. In all, we obtain a weighted average for to of 
1.61(9), which is also consistent with the value of 1.64 obtained 
by Engels et al. [15] in early study of Polyakov loop averages. 

For reference, we have included extrapolations of the crit- 
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2 We will present a more extensive treatment of the FSS of F tw in a forth- 
coming paper. 



Figure 5: Critical beta estimates vs N s relative to their infinite volume limits for 
N t = 4, 5 and 6, and the corresponding fits to (4). 
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Table 1 : Summary of results from the N% X N, lattices specified in columns 1 and 2 with the critical couplings for the largest two N s given explicitly in columns 3 
and 4; the infinite volume extrapolations from fits to Eq. (4) are shown in column 5 with the resulting exponents ai and^ 2 /dof. in columns 6 and 7. Columns 8 and 9 
show the same extrapolations when using the global average u> = 1.61(9) in all fits. Literature values are quoted for comparison from Refs. * * fl7ll and ' fljjll . 



ical coupling with u> fixed at 1.61. Due to correlations, the 
quoted errors for j3 c (N t )\ td= i,6\ should be taken with a grain of 
salt. Nevertheless, fixing u may lead to more accurate values 
of p c {N t ) for our large N, results. 

We furthermore obtain the fit parameter d{N,) in for each 
of the six N, values. This in turn allows us to fit its N, depen- 
dence. Using a two parameter form dfi t (N,) = d y Nj ca with a sin- 
gle effective exponent y& u a fairly good description is obtained 
with y flt = 3.97(7) (and d y = 0.06(1)). This fit works best for 
the smaller N t but it deteriorates somewhat towards N, -9. The 
N, dependence of d(N,) might well be determined by several 
competing terms with nearby exponents and is thus difficult to 
extract reliably from the data. Alternatively, using 3 parameter 
fits of the form 



d & (N,) = d s N? + d r Nj' 



(1) 



we obtain y^ t = 3.8(2) for 5 = 0, y^ t = 3.7(4) for 6 = co, or 
7flt = 3.5(6) for 8 = 0) + 1, for example, where we have used 
lo = 1.61 corresponding to our global fit for the correction to 
scaling exponent a>. Because <±> + 2 = 3.61, this suggests that 
one might also try a form 



dfc(N,) = d N? + di N? +1 + d 2 N? 



(8) 



with a) fixed at 1.61. This form is particularly interesting be- 
cause it will allow us below to parametrise the finite-size cor- 
rections entirely in terms of the aspect ratio in the form N t /N s . 
Unfortunately, the coefficients do and d\ of the two subleading 
powers in a> are much too correlated to determine all the 3 pa- 
rameters in this fit reliably from the available data. But we can 
get good fits with either d\ — or do = 0, see Fig. [6] Keeping 
the resulting two parameters fixed afterwards, a further fit to the 
remaining one (d\ or do) then reproduces a zero result with very 
high accuracy (i.e., values smaller than 10~ 4 in either case). 

It's interesting to check how much was gained by using 
the exact universal value (0 instead of pairwise intersections. 
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Figure 6: The parameter d(N,) from the fits of /3 C (N,, N s ) to Eq. (4) together 
with 2 parameter fits of its N, dependence via the form in Eq. (8) with d\ = 
(long-dashed) and do = (short-dashed), respectively. 



To this end we performed the pairwise intersection method for 
N t = 4. See Table [2] for the results. The left side of the ta- 
ble shows the intersection coupling for lattices with N s and 
N s = N s /2. On the right hand side are comparative results from 
intersections with the universal Ising reference line. 

It follows from Eq. Q that the pairwise intersection points 
should scale like Eq. (0J, except with -d — > 2(2" - l)d. As 
such, we've included an extrapolation with a> = 1.61 fixed. 

Since the simulations were not catered for pairwise intersec- 
tions, it may be unfair to directly compare the extrapolations 
in Table [2] Still, due to the rapid convergence of /3 c (N t ,N s ), 
our data was quite well centered around the intersection points 
except for the smallest lattices. And knowing what value of 
Ztw/Zg to concentrate the efforts around was a big advantage of 
our method. Anticipating the location of pairwise intersection 
points is much more troublesome. In all, it's clear that the ex- 
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Table 2: Comparison of critical couplings obtained for N, = 4 via pairwise 
intersection (left) and via intersection with the universal reference line (right). 
"This is the extrapolation with all 8 available N s values included, as quoted in 
the 8 th column of TablefT] The restriction to the 5 values of N s 6 (24, 96) listed 
here gives 6.53641(5), which is consistent but has an unnaturally small error. 



ploitation of the univer sal number F nv (/3 C ) = ln(l + 2 3/4 ) gave a 
significant boost to the precision of our results. 

4.2. 2D Ising model 

For the sake of comparison we repeated our procedure for 
the square 2D Ising model. Taking the exact solutions for Z ap 
and Z pp on NxN lattices from Ref. ifisill . we used Mathematica 
i2(ill to find 6 = ksT/J at the intersection of Z ap /Z pp with the 
universal value (0. Here J is the coupling of nearest neighbour 
spins. In Figure|7]we plot the results for N e [100,640]. The 
error bars are representative only (the errors are limited only by 
the working precision and should be smaller than 10~ 16 ). 

Again we fit the data with a FSS ansatz of the form 



em 



■cN 



-(w+l) 



In this case, we hold 6 C<<X fixed at the known value of 2/ ln(l + 
a/2). We obtain for the correction to scaling exponent u -2 + 6 
with 6 — » + as we increase the lower bound N m i n used in the 
fit, i.e., u> tends towards 2 from above {6 starts at around 2 • 10~ 4 
for the full range of N shown in Fig. [7] and it falls below 10~ 5 
at A^min around 400). Since the exponent of the leading irrel- 
evant operator that breaks rotational invariance is predicted to 
be exactly 2 112111 . our result is consistent with the conjecture 
of Ref. fl22ll that the only irrelevant operators that appear in the 
2D nearest neighbour Ising model are those due to the lattice 
breaking of rotational symmetry. This may be tested on a trian- 
gular 2D lattice where the leading irrelevant operator to break 
rotational invariance leads to a> — 4 instead, while the leading 
rotationally invariant operator would give an isotropic correc- 
tion to scaling with u — 2 in either case [22]. 

On the other hand, our correction exponent for SU(2) is 
clearly at odds with a> = 2. In this case, it's possible that there 
exists an irrelevant operator that is not present in the 2D Ising 
model or the corresponding conformal theory. It may be more 
likely, however, that our exponent is really an effective expo- 
nent. When there are several nearby competing exponents it 
is extremely difficult to extract the smallest one from simula- 
tions 
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Figure 7: Critical coupling estimate vs N for the 2D Ising model together with a 
fit by {9) yielding a> = 2.0002 for N mm = 100, or 01 = 2.00001 for JV min = 400. 



4.3. p c VS N, 

In 2+1 dimensions the coupling g 2 has the dimension of 
mass and sets the scale for the theory. The bare lattice coupling 
is then given by |fl7| 



2N C 



?3B = s\ + c\ag\ + c 2 a 2 g\ + 



Substituting this expansion into ft gives 

P 1 2^ 
= — j - ci - c 2 ag 3 + . . 
2N C ag\ 



(9) Note now that T = l/(N,a). So at criticality we have 



ftc(Nt) T Si 1 
= —r N> - C\ — Co H 

2N C g\ 2 T C N, 



(10) 



(11) 



(12) 



The first correction to the bare lattice coupling gives the y- 
intercept of /3 C (N,) vs N, while the second is a correction to 
linearity for small N,. 

In Fig. [8]we plot our results for ft c (N,), using the fitted val- 
ues with a) free. We include also the estimate yS e (3) — 4.978(35), 
which is the average of the literature values 4.943(13) iflill and 
5.013(15) ifTill . For the uncertainty we use their standard error, 
which is simply half their difference. The actual error used here 
has very little influence on the results. The data is fitted with 
both a straight line and with a function of the form ( fT2l . To 
the naked eye, it would appear that the linear approximation is 
good all the way down to N, = 3. However, our data is precise 
enough to pick up the deviation from linearity. In fact, the re- 
duced^- 2 of the linear fit is ~ 70. With, the 1 /N t correction the 



It was noted in 1151 that the observed correction to scaling exponent of 2+1- 



dimensional SU(2), at = 1.64 in their case or ai = 1.61(9) in ours, agreed well 
with some predictions for the universality class of the 2D Ising model which 
included 1.6 (23]. At the time it was discussed whether such non-integral cor- 
rection exponents could arise in other ferromagentic models of this class, and 
whether the corresponding correction amplitudes happened to vanish identi- 
cally in the exactly solvable pure Ising model, see [24] . This seems to be ruled 
out by a conformal field theory analysis: there is no irrelevant operator with 
oj < 2 in any unitary model of the 2D Ising class, see f25ll22Tl . 
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Figure 8: p c (Nt) vs JV,. Note the deviation from linearity. 



reduced^ 2 drops to 1.5. Thus, the data is very well described 
by Eq. ( fT2l >. We obtain the parameterisation 

Pc(N t )= 1.5028(21) JV,+ 0.705(21) - 0.718(49) i . (13) 

Consequently, from Eq. (fT2l for N c = 2, we obtain 



T c = 0.3757(5) g\. 



(14) 



This is compatible with the estimate T c /gi 0.385 with an 
error of +0.010 as quoted in Ref. uM- Note that it is the 1 /N, 
corrections included in (fT2l that lead to a somewhat lower value 
here as compared to 11811 which is however within their esti- 
mate of the systematic uncertainties. In order to translate our 
estimate into units of the zero-temperature string tension <x we 
use the mean value with standard error of the four N — 2 val- 
ues for y[crlg\ listed in B26I1 as pairs of upper and lower bounds 
including some systematic uncertainty, which gives ~\[<rlg\ = 
0.3347(5). With uncorrected error propagation, this together 
with our estimate (15[ then corresponds to 



TJyfcr = 1.1225(23), 



(15) 



which agrees very well with the corresponding result of lfl6ll . 
TJ^o- = 1.1224(90). 

From Eq. dot , the other constants in Eq. (fl2l) are analo- 
gously determined as 



-0.176(5), c 2 = 0.0675(5) . 



(16) 



Eq. ([TBI can be used to obtain accurate estimates of the critical 
coupling at large values of N,. We can furthermore include the 
finite-size corrections in our 'pseudo-critical coupling' (O for 
the intersection, at finite N s , of the vortex free energy F, H , in (0) 
with the universal value ln(l + 2 3 ^ 4 ). If we assume a form as 
given in Eq. (|S), these corrections can be expressed entirely in 
terms of the aspect ratio A = N,/N s . In particular, the leading 
corrections for small A can then conveniently be combined with 
the large N, expansion (fT2l as follows, 

p c (N t ,N s ) = (4^-d 2 A< J )N,-4c l -d l A< J (17) 

#3 

„2 , 

-Uc 2 — +d A p )— + ... , 

where we again used N c = 2 and p = o) + 1/v. Since v = 1 
in the 2D Ising model, our global fit to the correction to scaling 
exponent co amounts to p = co+\ = 2.61(9). If we assume 
d\ = 0, our fits yield 

d 2 = 0.1 34(4) , d = -0.5 1(9) , (18) 

and a reduced^- 2 of 2.7. With do — on the other hand, 

d 2 = 0.155(7) , d\ = -0.21(4) , (19) 

and a reduced^- 2 of 2.6. The corresponding fits are shown in 
Fig. [6] Apart from this systematic uncertainty in the other- 
wise elegant parametrisation © of the finite-size corrections, 
we have therefore determined all parameters in (TTTb . 
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4.4. f3 vs temperature 

So far, we have studied the critical couplings for different 
N,. These were measured from the intersection of the vortex 
free energy with the universal value. Since T = l/N t a, and 
with the temperature kept fixed at T c , this means that the corre- 
sponding lattice spacing at criticality scales inversely with N t , 
i.e., for two lattices with N, and N' t we have a' c /a c = N t /N' t , and 
at the given order in 1 /N, (in the infinite volume limit), 

^_^) = Z£ (iV ;_ i V t )- C2 ^(i 7 -i) + ....(20) 
4 4 g\ ' T c \N' t N,' 

Alternatively, we can consider a change of N, as a change in 
temperature at a fixed lattice spacing. In particular, a simulation 
done at criticality of the N' t lattice, with a c (N' t ) = a(N,), then 
corresponds to a simulation at T = {N',/N,)T C on the N, lattice, 
i.e., with /3 = /3(a), we have 

p c (N' t ) = p(\l(N' t T c )) = /3(l/(N,T)) . (21) 

Denoting this coupling by f3(T, N t ), we can therefore use Eq. ( f20b 
to write down an equation for the temperature dependence of 
the lattice coupling near criticality at a fixed N,, 

^_^) M £| ( 1_1 ) + ^ 

4 4 #3 N,^T T c > v ' 

If we compare simulations at different N, but with the same 
aspect ratio A = N,/N s in a finite volume, we can furthermore 
include the finite-size corrections here as in (fTTb - In terms of 
the reduced temperature t = T/T c - 1, we then have 

j3(T, N„ N s ) = p c (N„ N s ) + C(N„ N s ) t+ ... , (23) 

C(N„ N s ) = (4 ^ - d 2 A") N t + Uc 2 ^1+doAP)^-, 

V §3 ' V ic ' N, 

where fi c (N t , N s ) is the pseudo-critical coupling in Eq. (1171 1. and 
fl(T,N,,N s ) refers to the temperature dependence of the cou- 
pling near this /3 C at fixed N, and N s . The terms neglected are 
either subleading in the aspect ratio A, or of order 1 /N?, or they 
are of higher order in the reduced temperature. Note that criti- 
cality in a finite volume is a fuzzy concept, of course. As before, 
it is here defined by /3 = /3 C at t = from the intersection of the 
vortex free energy with the universal value, c.f., Eq. (0. 

Note that the linear relationship between lattice coupling ft 
and temperature in Eq. (l23l . which is of the form 

T/T c = l+(J3-/3 C )/C, (24) 

can be used to control the temperature at fixed N, and N s by 
adjusting the lattice coupling. The corresponding procedure for 
SU(2) in 3+1 dimensions, where /3 depends logarithmically on 
temperature, was used in H. We will exploit this to perform 
a detailed finite-size scaling analysis of the vortex free energy 
F tw in 2+1 dimensional SU(2) with high precision in a future 
paper. This will include scaling of data from various N, lattices. 



5. Summary 

We have shown that very accurate determinations of the 
critical couplings for the SU(2) deconfinement transition in 2+ 1 
dimensions are possible from intersecting vortex free energies 
with the known universal value from the 2D square Ising model 
at criticality. This allowed us to determine the critical coupling 
for lattices with up to N, -9 sites in the Euclidean time direc- 
tion. Its N, dependence could be determined with an accuracy 
sufficient to require significant 1 /N, corrections to the linear be- 
haviour of the continuum limit. As a result, the slope of this 
behaviour given by T c /gi might have been somewhat overesti- 
mated in the past. In particular, we observe that rather large N, 
may be required to reach the asymptotic continuum behaviour, 
even though it should be approached much more rapidly here 
than in 3+1 dimensions. 
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